Long-time tails in the random transverse Ising chain 
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Abstract 

Taking one-dimensional random transverse Ising model (RTIM) with the double- Gaussian disor- 
der for example, we investigated the spin autocorrelation function (SAF) and associated spectral 
density at high temperature by the recursion method. Based on the first twelve recurrants obtained 
analytically, we have found strong numerical evidence for the long-time tail in the SAF of a single 
spin. Numerical results indicate that when the standard deviation ajs (or ctbs) of the exchange 
couplings Jj (or the random transverse fields Bi) is small, no long-time tail appears in the SAF. 
The spin system undergoes a crossover from a central-peak behavior to a collective-mode behavior, 
which is the dynamical characteristics of RTIM with the bimodal disorder. However, when the 
standard deviation is large enough, the system exhibits similar dynamics behaviors to those of 
the RTIM with the Gaussian disorder, i.e., the system exhibits an enhanced central-peak behavior 
for large ajs or a disordered behavior for large (Jbs- In this instance, the long-time tails in the 
SAFs appear, i.e., C{t) ~ t^"^. Similar properties are obtained when the random variables (Jj or 
Bi) satisfy other distributions such as the double-exponential distribution and the double-uniform 
distribution. 
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The time-dependent behavior of quantum spin systems has been the subject of theoretical 
studies for quite some time. However, most hterature deals with the dynamic correlation 
functions of the pure quantum spin systems. It is found that the spin autocorrelation func- 
tions (SAFs) of these pure systems are associated with the temperature. The transverse 
correlation functions show a power-law behavior at T = OhJ, an exponential behavior at 

n nri 

< T < oo'2l, and a Gaussian behavior at T = oo [3|, |J] , respectively. As to the longitudinal 
correlation functions, Niemeijer 5| found that SAF of the XY spin chain is a damped oscil- 
lation function which may be expressed by the square of the Bessel function at T = oo, but 
it has a more complex expression for other temperature 6j. For the spin-1/2 van der Waals 
model, Dekeyser and Lee [7] pointed out that the power-law long-time tails are found in the 
transverse component only at high temperature. Only recently, since the disorder effects be- 
sides the effects of the temperature has been shown to affect the dynamics behavior of spin 
systems in a drastic way providing a very rich area of investigation, more and more atten- 
tion has been paid on the random quantum spin systems, especially on the low-dimensional 



systems, such as the random transverse Ising model (RTIM) js-ll^], XY chain ISHlTj] , XX 
and XXZ chain[l8|, etc.. The known results indicate that the disorder effects may drasti 



cally slow down the behavior of the disordered averaged spin-spin correlation function 18| , 
including the SAFs. For example, as the exchange couplings vary, the Gaussian behavior of 



xx-relaxation of the sing 



e impurity in XY chain slows down to the stretched exponential-like 



relaxation at T = cxo 19| . However, as far as we know, although the short-time behaviors 



of the SAFs have been reported 



10 



13 



15 



17l |. the long-time behaviors of them are still not 



understood. An important issue in this work is to explore whether it is possible to realize a 
situation in which the exchange couplings and the transverse fields are independent random 
variables, one can obtain the long-time behavior of SAF of the random transverse Ising 
chain. 

One-dimensional (1-D) RTIM is one of simplest but nontrivial example of the random 
quantum spin chains. It is defined by the following Hamiltonian: 



H 



-1/2 5^ J, 



1/2 5^ 5. 



(1) 



where 1/2 is the constant number introduced for the convenience of the derivation, af, af 
denote Pauli matrices at site i. The exchange couplings Jj (or the transverse fields Bi) are 
independent random variables obeying the distribution p(Jj) (or p{Bi)). This model has a 



well-known physical interpretation in connection not only with some quasi-one-dimensional 



hydrogen-bonded ferroelectric crystals like Cs(Hi xD.x)2P04, PbHi_^Da;P04 etc.[20|, |21| . 



lis model may capture vital 



14 



23| 



but also with Ising spin glasses LiHoo.i67Yo.833F4[22|. Besides, t 
features of the recent neutron scattering experiment in LiHoF4 

In this rapid communication, we investigate the dynamics of RTIM with the double- 
Gaussian disorder, i.e., Ji or Bi respectively satisfy the double-Gaussian distribution 



p(/3.)=m(/3i) + (i-p)P2(A 



(2) 



where 



and 



Pi(A) 



1 



27ro"/ 



exp [-(/3^-/3^jV2aJ] 



exp [-(/3^-/32jV2a^] 



(3) 



(4) 



are the standard Gaussian distributions, wherein Pi^^ {Jim or Bim) and /32m {-hm or B2m) 
denote the corresponding mean values of the random variables /3j (Jj or {aj or 

as) is the standard deviation. The physical sense of Eq. ([2]) is that the random variable 
/3j satisfies Pi(/3j) with probability p and P2(/^j) with probability (1 — p), respectively. The 
mean value and the standard deviation of the random variables (3^ obeying distribution p(/3j) 
are respectively p/J^^ + {l-p)(32m and a/^s = ^^p{l - p){,Pim - P2mY + Since a^s is an 
increasing function of o"^, for convenience, ap is used to denote a^s in the following unless 
specified otherwise. 



Obviously, the cases discussed in references [10| and are two special cases of the 
RTIM with the double- Gaussian disorder. When ap — )• 0, Eq. ([2]) becomes the bimodal 
distribution 

p(/3J = p - + {l-p) m - P2m)- (5) 

When p = 1 (or p = 0), Eq. ([2]) becomes a standard Gaussian distribution. Of course, Pi(/3j) 
and P2{Pi) may be other distributions, such as an exponential distribution or a uniform 
distribution. In this way, p(/3j) becomes the double-exponential distribution or the double- 
uniform distribution. 

The basic tool employing in this work is the recursion method 
in the study on many-body dynamics P 
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11 



15 



17 



25 



24( 1 , which is very powerful 



29|. It may help us obtain the 



3 



SAF and associated spectral density. They are respectively defined as 



C{t) = {a^{t)a^{0)), (6) 

and 

r+oo 

$ (w) = / dtC (t) e-^'^*, (7) 







where (■ ■ ■ ) denotes an average over the random variables is performed after the statistical 
average. At T — )■ oo, the SAF may be expanded into a power series 

where is the 2A;th moment of the SAF. And the spectral density may be determined via 
the relation 

$ (u) = lim Reao{z) \z=e+iu, (9) 



e-i-O 



where 



ao{z) = l/{z + Ai/ {z + A2/ [z + A3/ {z + ...)]}), (10) 

which is the Laplace transform of the SAF ([6]). The continued- fraction coefficients (recur- 
rants) ^ui = 1, 2, ■ ■ ■ ) of the SAF may be derived by the recurrence relation I given in 



reference 24] • 



Although only a finite number of recurrants can be deduced analytically, they play vital 
roles in the whole study process. On the one hand, since the first 2k moments may be 
expressed by the first k recurrants [24], SAF can be determined by the Fade approximant 
constructed out of the known moments. On the other hand, one can determine the spectral 

n 



density 



24 



28 



of the SAF by employing the known recurrants and so-called Gaussian terminator |lOl. 



30| . Due to the first 12 recurrants have been deduced analytically, we extract more 



ll| in which 



relevant information from them in this paper than we did in our previous work 
only the first 9 recurrants have been obtained analytically. 

In this work we analyze two types of RTIM: a random-bonds model, and a random-fields 
model. In the former, the exchange couplings Jj are chosen independently from the double- 
Gaussian distributions defined by Eq. (jj]), while the transverse fields Bj remain unaltered. 



To compare our results with the previous results presented in reference lOj, we set the same 
parameters as those in it: Jim = 1.0, = 0.4 and Bi = B = 1 which fixes the energy scale, 
and a J varies from 0.01 to 2. C {t) and $ {u) are plotted in Fig. 1 and Fig. 2, respectively. 



Each inset in Fig. 2 gives the recurrants A,y (z/ = 1, ■ ■ ■ , 12) for the same parameters in each 
figure, and the fines in it are just a guide to the eye. 

As a J — )■ 0, the double- Gaussian distribution turns into the bimodal distribution. Taking 
a J = 0.01 for example, as shown in Figs. 1 (a) and 2 (a), the dynamics of the system under- 
goes a crossover from a central-peak behavior to a collective-mode behavior as p decreases 
from 1 to 0. Although we can not obtain C (t) for large time, comparison shows that Fig.l 
(a) offers a better improvement in SAFs for short time than those given by Florencio, et 



al. 



lOj. Comparing all the figures in Fig. 1 and Fig. 2, one finds that as crj increases from 



0.01 to 2, the damped oscillations of the SAFs diminish in magnitude, meanwhile the recur- 
rants for the same parameters tend to fall into a straight line. In this process, the spectral 
weights at w = are enhanced, and the peaks' locations of the spectral lines move towards 
u = gradually (e.g., line p = 0). The solid lines with p = 1 in Fig. 2 illustrate how an 
almost pure Gaussian profile of ^{u) evolves into a curve with some additional structure 
which consists of a central peak of increasing height and decreasing width. Above analysis 
indicates that the collective-mode behavior becomes weaker and weaker, while the dynamics 
becomes increasingly dominated by a central-peak behavior as aj increases. When aj are 
large enough (see Figs. 1(d) and 2(d)), the crossover vanishes, and the system exhibits an 
enhanced central-peak behavior. These results are in agreement with those in our previous 



work 



The infrared (i.e., low-frequency) singularity in the spectral density signals the SAF 
should manifest itself a characteristic quantum-diffusion long-time tailSlI- We find, as 
expected, long-time tails exist in the SAFs when aj are large enough. Fig. 3 gives the log- 
log plot of the absolute values of C (t) versus time. The solid lines with p = 1 in Figs. 1 and 
3 show how a Gaussian behavior of C {t), due to the disorder effects, slows down to a power- 
law relaxation. Calculations indicate that C (t) ~ for t > 20. This result is identical 
with the pure XY interaction case {R = 0) of the spin-1/2 van der Waals model[3]. In fact, 
the XY model and the transverse Ising model are dynamically equivalent to the spin van 
der Waals modelj4|. These systems are indicated by the same geometry of their respective 
dynamics Hilbert spaces. Finally, we should like to emphasize that success in obtaining the 
long-time tails of the SAFs is owing to the fact that the true long-time behavior is only 



nebulously encoded in the first few continued-fraction coefficients 



3l[. 



For the random-fields model, the random transverse fields Bi are independent and satisfy 
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the double-Gaussian distributions given by Eq. ([2]), while the exchange couplings are uni- 
form (i.e., Ji = J = 1). For convenience to compare our results with the known results 10[ |. 
we set Bim = 0.6, = 1-4, and as varies from 0.01 to 2. The SAFs and the correspond- 
ing spectral densities are shown in Fig. 4 and Fig. 5, respectively. The inserts in them 
respectively give the double-logarithmic plots of the absolute values of C (t) versus time and 
the associated recurrants (z/ = 1, ■ ■ ■ , 12) for the same parameters in each figure. 

When (Jb is small (e.g., cr^ < 0.10), the system undergoes a crossover from a central-peak 
behavior to a collective-mode behavior as p decreases from 1 to 0, as shown in Figs. 4 (a) 
and 5 (a). This result is consistent with that obtained by Florencio, et al. 10|. However, as 
as increases, this crossover disappears gradually, see Fig. 5 (b), (c) and (d). Meanwhile, 
the SAFs exhibit an interesting phenomenon. From Fig. 4 (b) we find that the decays of 
long-time tails in lines p = and p = 1 are proportional to — and those in other lines 
are proportional to t"^. When is large enough (e.g., ctb = 1.00), SAFs fall into three 
successive phases. In the initial phase (0 < t < 1), the decay of C (t) is Gaussian regime. In 
the third phase [t > 20), the long-time tails of SAFs obey the power- law decay There is 
a cross-over time between the Gaussian regime and the power-law regime. In this cross-over, 
the oscillation of C {t) implies that the x component of the spin may be coherently coupled 
to some mode with small fluctuation [3]. Thus, slow decay may emerge. The strong peaks 
at a; = aslo signal the existence of the long-time tails in the SAFs. The cross-over time 
becomes shorter as cr^ is increased further. While more and more high-frequency modes are 
involved in the corresponding spectral densities. The system is in a disordered state when 
(Tb is large enough. 

Above we have investigated the dynamics of the RTIM with the double-Gaussian disor- 
der at high temperature by the recursion method. This generalized model may recover the 
dynamics of both the RTIM with the bimodal disorder and that with the Gaussian disorder. 
One may wonder if these spin dynamics features are universal or not. Therefor, we have 
studied the dynamics of the RTIM with other types of disorder, in which the variables inde- 
pendently satisfy the double-Exponential distribution, and the double-Uniform distribution. 
Numerical results indicate that these systems show a similar dynamical behavior to those 
of the RTIM with the double-Gaussian disorder. Generally speaking, when ajs (or cbs) is 
small, the dynamics of the system undergoes a crossover from a central-peak behavior to a 
collective-mode one as p (or the mean value of Ji or Bi ) varies. No long-time tail in SAF 



6 



is obtained. As ajs (or ubs) increases, the crossover vanishes gradually, and the power- law 
decay appears in the SAFs when the standard deviation is large enough. For the case of 
large ajs the system exhibits a central-peak behavior, but it shows a disordered behavior 
for the case of large (Jbs- As expected, we found that, independently of the type of disorder, 
the SAFs in 1-D RTIM decay as for large time due to quantum spin diffusion. 
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Figure Captions 



Fig. 1. (Color online) Spin autocorrelation functions C (t) of RTIM with the double- 
Gaussian random bonds. Better improvement in SAFs is made than those given in references 







m and 



Fig. 2. (Color online) The associated spectral densities for the same parameters as in 
Fig. 1. The insets show the recurrants A^, {u = 1, ■ ■ ■ ,12). The lines in the insets are 
just a guide to the eye. For small aj (e.g., aj < 0.01), (a) recovers $(a;) of RTIM with 
the bimodal distribution. As sigmaj increases, the crossover from a central-peak behavior 
to a collective-mode behavior vanishes gradually. The system only exhibits an enhanced 
central-peak behavior when aj is large enough. 

Fig. 3. (Color online) Log-log plot of \C (t)\ versus t for different values of p and aj. 
Black dashed line in each figure represents a power-law decay t~^. 

Fig. 4. (Color online) SAFs C (t) of RTIM with the double- Gaussian random fields. 
When aB is small (e.g., = 0.01), (a) approximately recovers the results presented in Ref. 
[lO|. The inserts present log-log plot of \C {t)\ versus t. The black dashed line in each insert 
represents a power-law decay t~^. 

Fig. 5. (Color online) The corresponding spectral densities for the same parameters as in 
Fig. 4. For small aB (e.g., ctb = 0.01), the system moves from a collective mode dynamics 
to a central peak type of dynamics as p varies from to 1. However, the system is in a 
disordered state when is large enough (e.g., = 2.00). 
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